3.1.52 \(\int \frac {1}{(d+e x^n)^2 (a+c x^{2 n})^2} \, dx\) [52]

3.1.52.1 Optimal result
3.1.52.2 Mathematica [A] (verified)
3.1.52.3 Rubi [A] (verified)
3.1.52.4 Maple [F]
3.1.52.5 Fricas [F]
3.1.52.6 Sympy [F(-1)]
3.1.52.7 Maxima [F]
3.1.52.8 Giac [F]
3.1.52.9 Mupad [F(-1)]

3.1.52.1 Optimal result

Integrand size = 21, antiderivative size = 410 \[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\frac {c x \left (c d^2-a e^2-2 c d e x^n\right )}{2 a \left (c d^2+a e^2\right )^2 n \left (a+c x^{2 n}\right )}+\frac {c e^2 \left (3 c d^2-a e^2\right ) x \operatorname {Hypergeometric2F1}\left (1,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a \left (c d^2+a e^2\right )^3}-\frac {c \left (c d^2-a e^2\right ) (1-2 n) x \operatorname {Hypergeometric2F1}\left (1,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{2 a^2 \left (c d^2+a e^2\right )^2 n}+\frac {4 c e^4 x \operatorname {Hypergeometric2F1}\left (1,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )}{\left (c d^2+a e^2\right )^3}-\frac {4 c^2 d e^3 x^{1+n} \operatorname {Hypergeometric2F1}\left (1,\frac {1+n}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a \left (c d^2+a e^2\right )^3 (1+n)}+\frac {c^2 d e (1-n) x^{1+n} \operatorname {Hypergeometric2F1}\left (1,\frac {1+n}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a^2 \left (c d^2+a e^2\right )^2 n (1+n)}+\frac {e^4 x \operatorname {Hypergeometric2F1}\left (2,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )}{d^2 \left (c d^2+a e^2\right )^2} \]

output
1/2*c*x*(c*d^2-a*e^2-2*c*d*e*x^n)/a/(a*e^2+c*d^2)^2/n/(a+c*x^(2*n))+c*e^2* 
(-a*e^2+3*c*d^2)*x*hypergeom([1, 1/2/n],[1+1/2/n],-c*x^(2*n)/a)/a/(a*e^2+c 
*d^2)^3-1/2*c*(-a*e^2+c*d^2)*(1-2*n)*x*hypergeom([1, 1/2/n],[1+1/2/n],-c*x 
^(2*n)/a)/a^2/(a*e^2+c*d^2)^2/n+4*c*e^4*x*hypergeom([1, 1/n],[1+1/n],-e*x^ 
n/d)/(a*e^2+c*d^2)^3-4*c^2*d*e^3*x^(1+n)*hypergeom([1, 1/2*(1+n)/n],[3/2+1 
/2/n],-c*x^(2*n)/a)/a/(a*e^2+c*d^2)^3/(1+n)+c^2*d*e*(1-n)*x^(1+n)*hypergeo 
m([1, 1/2*(1+n)/n],[3/2+1/2/n],-c*x^(2*n)/a)/a^2/(a*e^2+c*d^2)^2/n/(1+n)+e 
^4*x*hypergeom([2, 1/n],[1+1/n],-e*x^n/d)/d^2/(a*e^2+c*d^2)^2
 
3.1.52.2 Mathematica [A] (verified)

Time = 0.47 (sec) , antiderivative size = 298, normalized size of antiderivative = 0.73 \[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\frac {x \left (\frac {c e^2 \left (3 c d^2-a e^2\right ) \operatorname {Hypergeometric2F1}\left (1,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a}+4 c e^4 \operatorname {Hypergeometric2F1}\left (1,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )-\frac {4 c^2 d e^3 x^n \operatorname {Hypergeometric2F1}\left (1,\frac {1+n}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a (1+n)}+\frac {c \left (c d^2-a e^2\right ) \left (c d^2+a e^2\right ) \operatorname {Hypergeometric2F1}\left (2,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a^2}+\frac {e^4 \left (c d^2+a e^2\right ) \operatorname {Hypergeometric2F1}\left (2,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )}{d^2}-\frac {2 c^2 d e \left (c d^2+a e^2\right ) x^n \operatorname {Hypergeometric2F1}\left (2,\frac {1+n}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a^2 (1+n)}\right )}{\left (c d^2+a e^2\right )^3} \]

input
Integrate[1/((d + e*x^n)^2*(a + c*x^(2*n))^2),x]
 
output
(x*((c*e^2*(3*c*d^2 - a*e^2)*Hypergeometric2F1[1, 1/(2*n), (2 + n^(-1))/2, 
 -((c*x^(2*n))/a)])/a + 4*c*e^4*Hypergeometric2F1[1, n^(-1), 1 + n^(-1), - 
((e*x^n)/d)] - (4*c^2*d*e^3*x^n*Hypergeometric2F1[1, (1 + n)/(2*n), (3 + n 
^(-1))/2, -((c*x^(2*n))/a)])/(a*(1 + n)) + (c*(c*d^2 - a*e^2)*(c*d^2 + a*e 
^2)*Hypergeometric2F1[2, 1/(2*n), (2 + n^(-1))/2, -((c*x^(2*n))/a)])/a^2 + 
 (e^4*(c*d^2 + a*e^2)*Hypergeometric2F1[2, n^(-1), 1 + n^(-1), -((e*x^n)/d 
)])/d^2 - (2*c^2*d*e*(c*d^2 + a*e^2)*x^n*Hypergeometric2F1[2, (1 + n)/(2*n 
), (3 + n^(-1))/2, -((c*x^(2*n))/a)])/(a^2*(1 + n))))/(c*d^2 + a*e^2)^3
 
3.1.52.3 Rubi [A] (verified)

Time = 0.63 (sec) , antiderivative size = 410, normalized size of antiderivative = 1.00, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.095, Rules used = {1767, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {1}{\left (a+c x^{2 n}\right )^2 \left (d+e x^n\right )^2} \, dx\)

\(\Big \downarrow \) 1767

\(\displaystyle \int \left (-\frac {c e^2 \left (a e^2-3 c d^2+4 c d e x^n\right )}{\left (a e^2+c d^2\right )^3 \left (a+c x^{2 n}\right )}-\frac {c \left (a e^2-c d^2+2 c d e x^n\right )}{\left (a e^2+c d^2\right )^2 \left (a+c x^{2 n}\right )^2}+\frac {4 c d e^4}{\left (a e^2+c d^2\right )^3 \left (d+e x^n\right )}+\frac {e^4}{\left (a e^2+c d^2\right )^2 \left (d+e x^n\right )^2}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {c^2 d e (1-n) x^{n+1} \operatorname {Hypergeometric2F1}\left (1,\frac {n+1}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a^2 n (n+1) \left (a e^2+c d^2\right )^2}-\frac {c (1-2 n) x \left (c d^2-a e^2\right ) \operatorname {Hypergeometric2F1}\left (1,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{2 a^2 n \left (a e^2+c d^2\right )^2}-\frac {4 c^2 d e^3 x^{n+1} \operatorname {Hypergeometric2F1}\left (1,\frac {n+1}{2 n},\frac {1}{2} \left (3+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a (n+1) \left (a e^2+c d^2\right )^3}+\frac {c e^2 x \left (3 c d^2-a e^2\right ) \operatorname {Hypergeometric2F1}\left (1,\frac {1}{2 n},\frac {1}{2} \left (2+\frac {1}{n}\right ),-\frac {c x^{2 n}}{a}\right )}{a \left (a e^2+c d^2\right )^3}+\frac {c x \left (-a e^2+c d^2-2 c d e x^n\right )}{2 a n \left (a e^2+c d^2\right )^2 \left (a+c x^{2 n}\right )}+\frac {4 c e^4 x \operatorname {Hypergeometric2F1}\left (1,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )}{\left (a e^2+c d^2\right )^3}+\frac {e^4 x \operatorname {Hypergeometric2F1}\left (2,\frac {1}{n},1+\frac {1}{n},-\frac {e x^n}{d}\right )}{d^2 \left (a e^2+c d^2\right )^2}\)

input
Int[1/((d + e*x^n)^2*(a + c*x^(2*n))^2),x]
 
output
(c*x*(c*d^2 - a*e^2 - 2*c*d*e*x^n))/(2*a*(c*d^2 + a*e^2)^2*n*(a + c*x^(2*n 
))) + (c*e^2*(3*c*d^2 - a*e^2)*x*Hypergeometric2F1[1, 1/(2*n), (2 + n^(-1) 
)/2, -((c*x^(2*n))/a)])/(a*(c*d^2 + a*e^2)^3) - (c*(c*d^2 - a*e^2)*(1 - 2* 
n)*x*Hypergeometric2F1[1, 1/(2*n), (2 + n^(-1))/2, -((c*x^(2*n))/a)])/(2*a 
^2*(c*d^2 + a*e^2)^2*n) + (4*c*e^4*x*Hypergeometric2F1[1, n^(-1), 1 + n^(- 
1), -((e*x^n)/d)])/(c*d^2 + a*e^2)^3 - (4*c^2*d*e^3*x^(1 + n)*Hypergeometr 
ic2F1[1, (1 + n)/(2*n), (3 + n^(-1))/2, -((c*x^(2*n))/a)])/(a*(c*d^2 + a*e 
^2)^3*(1 + n)) + (c^2*d*e*(1 - n)*x^(1 + n)*Hypergeometric2F1[1, (1 + n)/( 
2*n), (3 + n^(-1))/2, -((c*x^(2*n))/a)])/(a^2*(c*d^2 + a*e^2)^2*n*(1 + n)) 
 + (e^4*x*Hypergeometric2F1[2, n^(-1), 1 + n^(-1), -((e*x^n)/d)])/(d^2*(c* 
d^2 + a*e^2)^2)
 

3.1.52.3.1 Defintions of rubi rules used

rule 1767
Int[((d_) + (e_.)*(x_)^(n_))^(q_)*((a_) + (c_.)*(x_)^(n2_))^(p_), x_Symbol] 
 :> Int[ExpandIntegrand[(d + e*x^n)^q*(a + c*x^(2*n))^p, x], x] /; FreeQ[{a 
, c, d, e, n, p, q}, x] && EqQ[n2, 2*n] && NeQ[c*d^2 + a*e^2, 0] && ((Integ 
ersQ[p, q] &&  !IntegerQ[n]) || IGtQ[p, 0] || (IGtQ[q, 0] &&  !IntegerQ[n]) 
)
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.52.4 Maple [F]

\[\int \frac {1}{\left (d +e \,x^{n}\right )^{2} \left (a +c \,x^{2 n}\right )^{2}}d x\]

input
int(1/(d+e*x^n)^2/(a+c*x^(2*n))^2,x)
 
output
int(1/(d+e*x^n)^2/(a+c*x^(2*n))^2,x)
 
3.1.52.5 Fricas [F]

\[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\int { \frac {1}{{\left (c x^{2 \, n} + a\right )}^{2} {\left (e x^{n} + d\right )}^{2}} \,d x } \]

input
integrate(1/(d+e*x^n)^2/(a+c*x^(2*n))^2,x, algorithm="fricas")
 
output
integral(1/(a^2*e^2*x^(2*n) + 2*a^2*d*e*x^n + a^2*d^2 + (c^2*e^2*x^(2*n) + 
 2*c^2*d*e*x^n + c^2*d^2)*x^(4*n) + 2*(a*c*e^2*x^(2*n) + 2*a*c*d*e*x^n + a 
*c*d^2)*x^(2*n)), x)
 
3.1.52.6 Sympy [F(-1)]

Timed out. \[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\text {Timed out} \]

input
integrate(1/(d+e*x**n)**2/(a+c*x**(2*n))**2,x)
 
output
Timed out
 
3.1.52.7 Maxima [F]

\[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\int { \frac {1}{{\left (c x^{2 \, n} + a\right )}^{2} {\left (e x^{n} + d\right )}^{2}} \,d x } \]

input
integrate(1/(d+e*x^n)^2/(a+c*x^(2*n))^2,x, algorithm="maxima")
 
output
(c*d^2*e^4*(5*n - 1) + a*e^6*(n - 1))*integrate(1/(c^3*d^8*n + 3*a*c^2*d^6 
*e^2*n + 3*a^2*c*d^4*e^4*n + a^3*d^2*e^6*n + (c^3*d^7*e*n + 3*a*c^2*d^5*e^ 
3*n + 3*a^2*c*d^3*e^5*n + a^3*d*e^7*n)*x^n), x) - 1/2*(2*(c^2*d^2*e^2 - a* 
c*e^4)*x*x^(2*n) + (c^2*d^3*e + a*c*d*e^3)*x*x^n - (c^2*d^4 - a*c*d^2*e^2 
+ 2*a^2*e^4)*x)/(a^2*c^2*d^6*n + 2*a^3*c*d^4*e^2*n + a^4*d^2*e^4*n + (a*c^ 
3*d^5*e*n + 2*a^2*c^2*d^3*e^3*n + a^3*c*d*e^5*n)*x^(3*n) + (a*c^3*d^6*n + 
2*a^2*c^2*d^4*e^2*n + a^3*c*d^2*e^4*n)*x^(2*n) + (a^2*c^2*d^5*e*n + 2*a^3* 
c*d^3*e^3*n + a^4*d*e^5*n)*x^n) - integrate(1/2*(a^2*c*e^4*(4*n - 1) - c^3 
*d^4*(2*n - 1) - 6*a*c^2*d^2*e^2*n + 2*(a*c^2*d*e^3*(5*n - 1) + c^3*d^3*e* 
(n - 1))*x^n)/(a^2*c^3*d^6*n + 3*a^3*c^2*d^4*e^2*n + 3*a^4*c*d^2*e^4*n + a 
^5*e^6*n + (a*c^4*d^6*n + 3*a^2*c^3*d^4*e^2*n + 3*a^3*c^2*d^2*e^4*n + a^4* 
c*e^6*n)*x^(2*n)), x)
 
3.1.52.8 Giac [F]

\[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\int { \frac {1}{{\left (c x^{2 \, n} + a\right )}^{2} {\left (e x^{n} + d\right )}^{2}} \,d x } \]

input
integrate(1/(d+e*x^n)^2/(a+c*x^(2*n))^2,x, algorithm="giac")
 
output
integrate(1/((c*x^(2*n) + a)^2*(e*x^n + d)^2), x)
 
3.1.52.9 Mupad [F(-1)]

Timed out. \[ \int \frac {1}{\left (d+e x^n\right )^2 \left (a+c x^{2 n}\right )^2} \, dx=\int \frac {1}{{\left (a+c\,x^{2\,n}\right )}^2\,{\left (d+e\,x^n\right )}^2} \,d x \]

input
int(1/((a + c*x^(2*n))^2*(d + e*x^n)^2),x)
 
output
int(1/((a + c*x^(2*n))^2*(d + e*x^n)^2), x)